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, We have shown that the e-entropy determined from the time-series of the velocity fluctuation in 

f"*^ ' a forced turbulence simulated by the lattice Boltzmann method obeys the Kolmogorov scaling. 
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Turbulence is one of the typical phenomena with multiscalc motions. Each phenomena of a scale strongly couples 
with all the other scales of turbulent motion. In order to analyze such a system a scale-dependent entropy, the so-called 
e-entropy, works well. Q For example, the time series of the velocity fluctuation in turbulence leads to a non-trivial 
i—) | scaling relation of the e-entropy, h(e) <x e -3 , expected from the Kolmogorov scaling. (2) The existing experimental 
data are consistent with this scaling. p| In this Short Note we try to show the consistency of the scaling with the 
simulation on the basis of the lattice Boltzmann method which is one of the easiest ways to simulate turbulence. 

The distribution function f%(r,t) of the particle with the velocity q at position r and time t obeys the lattice 
Boltzmann equation MM 



fi(r + 4 t + 1) - Mr, t) = -u[f t (r, t) - /°(f, t)}, (1) 

r/3 ; 

I where we have adopted the Bhatnagar-Gross-Krook (BGK) model for the collision term and the time step has been 
i chosen as unity. The position r of the particle is restricted on the cubic lattice and the lattice spacing is chosen as 
unity. The local equilibrium distribution ff(r,t) is assumed to be given as 

I3 i f°(r, t) = w lP [l + 3@ ■ u) + §@ • uf - § (u ■ u)}. (2) 
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The density p and the velocity u of the fluid are given by the sum of the contributions of the particles as p = fi{r, t) 
and pu = fi(r,t)ci, respectively. Using the multiscale analysis the Navier-Stokes equation, which is supposed to 
be able to describe turbulence, is derived from the present BGK lattice Boltzmann equation and the corresponding 
viscosity v — (l/u> — l/2)/3. 

In the following we use a 15- velocity model and is chosen as C\ = (0, 0, 0), c j; = (±1, 0, 0), (0, ±1, 0), (0, 0, ±1) for 
% = 2, 3, 7 and Cj = (±1, ±1, ±1) for « = 8, 9, — , 15. Then the weight factor Wi for the local equilibrium distribution 
is determined as w\ = 2/9, Wi — 1/9 for i = 2, 3, ■ • •, 7 and Wi = 1/72 for i = 8, 9, • • •, 15. The relaxation frequency u 
■ is chosen asw = 1.94. 

{SJ ' In order to simulate forced turbulence we add the forcing term <fj • F(r, t) due to an applied body force F(r, t) to 
the right hand side of eq. (1). In order to reproduce the Navier-Stokes equation we should choose as gt — ci/10 for 
the present 15- velocity model. Q For simplicity we apply a solenoidal force in y-direction, F = (0, F y (z, t),0), and 
F y (z, t) is a Gaussian white noise whose variance a is a function of the lattice coordinate in z-direction, a = |0.01 x 
sin(27r/c/2;/L)|, where L is the linear dimension of the simulation region of cubic box and L = 100 (1 < x,y, z < 100). 
^ , We have adopted the periodic boundary condition. Although we have chosen as fc/ = 4, the results of the simulation 
£H ■ is insensitive to the choice. 

In the initial state of the simulation the fluid density is uniform, p = 1, and the spatial distribution of the fluid 
O ■ velocity u is Gaussian with the variance a = V0.0075. 



The Reynolds number R for our simulation [g|4J| is estimated as R = y/ (u 2 )L/2-kv and time-dependent as shown 
, in Fig. 1. Here (• • •) represents the spatial average. In our simulation we can realize a turbulent state for relatively 
/\ " small Reynolds number, since we add random force. 

We focus our attention to the time series of the y-component u y (r,t) of the fluid velocity u(r,t) as shown in Fig. 
2, since our system is anisotropic due to the forcing. This turbulent signal is characterized by the exit-time r(e; r) at 
position r. At the exit-time t = r(e;r) the first exit satisfying the condition, \u y (r, to + i) — u y (r,to)\ > e/2, occurs. 
The time t is measured from to and to = 1000 in our simulation regarding the states at the first 1000 time-steps as 
transitional. 

The mean of the exit-time, r(e) = (r(e; r)), is related to the e-entropy, h(e), as Q 

K e ) K ^y- ( 3 ) 
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In the region where Kolmogorov's 4/5 law holds g|| the velocity difference between two points behaves as (\u(f + 
R,to) — u(r, to)]) oc l^l 1 / 3 . Under the assumption of Taylor's hypothesis, H|| [u(r, t + 1) — u(r, £ )] ~ [u(f+R,t )~ 
u{r, tg)], with ]R] = Ut, where U is the mean velocity for large scale sweeping, we obtain 

(\u(r,t +t)-u(r,t )\) «( 1/3 , (4) 

as a mean- field description. Thus the definition of the exit-time leads to r(e) 1 / 3 oc e. From eq. (3) we can conclude 
that h(e) oc e~ 3 in the above mean-field sense. 

In accordance with the above scaling analysis our simulation data for l/r(e) shown in Fig. 3 have the scaling region 
for e > eo where eo ~ 0.04. Thus by simulation we have established that h(e) oc e~ 3 . For e < eo, namely for short-time 
events, we can not expect the scaling, since the mean-field assumption for large scale sweeping is valid for the average 
of many events and not applicable to a few events at short time- interval. These behaviors are consistent with the 
experimantal data. For e< eo the exit in the simulation occurs within the time step. 

In summary we have shown that the velocity fluctuation in a forced turbulence simulated by the lattice Boltzmann 
method leads to the e-entropy scaling, h(e) oc e~ 3 , expected from the Kolmogorov scaling. Although the scaling 
relation is a robust result of the mean-field description, the intermittent aspects of turbulence are beyond the scope 
of the present study. 
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FIG. 1. The Reynolds number as a function of time. 

FIG. 2. The y-component of the fluid velocity u y as a function of time at r — (50, 50, 50). 

FIG. 3. The inverse of the mean exit-time r(e). Here the dots and straight line represent the simulation data and e -3 law, 
respectively. 



2 




500 1000 1500 2000 2500 3000 



t 



